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1^-^ , Abstract. Based on recent work on spectral decomposition of the emission of star-forming galaxies, we assess whether 

■ the integrated 2-10 keV emission from high-mass X-ray binaries (HMXBs), L^^{q^, can be used as a reliable estimator of 

' ongoing star formation rate (SFR). Using a sample of 46 local (z^ 0.1) star-forming galaxies, and spectral modeling of ASCA, 

, BeppoSAX, and XMM-Newton data, we demonstrate the existence of a linear SFR-L^_^^^ relation which holds over ~5 

' decades in X-ray luminosity and SFR. The total 2-10 keV luminosity is not a precise SFR indicator because at low SFR (i.e., 

I in normal and moderately-starbursting galaxies) it is substantially affected by the emission of low-mass X-ray binaries, which 

do not trace the current SFR due to their long evolution lifetimes, while at very high SFR (i.e., for very luminous FIR-selected 
I ' galaxies) it is frequently affected by the presence of strongly obscured AGNs. The availability of purely SB-powered galaxies 

O I - whose 2-10 keV emission is mainly due to HMXBs - allows us to properly calibrate the SFR-L^.^'J^'^ relation. The SFR- 

L^^io^ relation holds also for distant (z ~ 1) galaxies in the Hubble Deep Field North sample, for which we lack spectral 
^ ' information, but whose SFR can be estimated from deep radio data. If confirmed by more detailed observations, it may be 

possible to use the deduced relation to identify distant galaxies that are X-ray overluminous for their (independently estimated) 
SFR, and are therefore likely to hide strongly absorbed AGNs. 

X 
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1 . Introduction these classes of source. The spectral properties and relative 

abundances of the various classes of stellar sources determine 
X-ray emission of star-formins galaxies (SFGs) consists of ^, -^r , • r ^ n w 

■' o o V / jjjg composite X-ray spectrum arising from a stellar population 

various components including discrete sources, such as X- ^ , rr. ^ nnrv-i " . n " 

... . ,„.,„ ,..„„ . of Galactic composition. Therefore, the PR02 stellar contn- 



ray binaries and supernova remnants (SNRs), diffuse hot gas 
Compton scattering of ambient FIR photons, and possibly an 
active nucleus. The resulting integrated spectra harbor the sig- 



bution has then no essential degrees of freedom: fixed by the 

synthetic model of the Galactic population of binaries and by 

. . the observed X-ray spectra of the contributing components, it 

natures of these emission components. . v . •„ j i. t ■ c 

^ represents the X-ray spectrum emitted by a Galactic mix of 

Persic & Rephaeli (2002, hereafter PR02) have quan- j^^^Bs and LMXBs. As such, the stellar component of the 

titatively assessed the roles of the various X-ray emission nnm . i . . • i i i . u ■ ^ c 

■' ■' PR02 template spectrum is likely to be appropriate for a qui- 

mechanisms in SFGs. They have used an equilibrium stellar- -.i ^ • i i i -.u \/r n wr i i i 

■' ^ etly star forming galaxy like the Milky Way (see Fig.l-/e/f). 

population synthesis model of the Galactic population of X-ray i 

binaries (Iben et al. 1995a,b) to deduce birthrates for interact- 
ing binaries; these, combined with estimates of the duration of 

the X-ray bright phase, have allowed PR02 to make reaUstic es- ' The diffuse non-stellar pai't of the PR02 synthetic spectrum has 

timates of the relative (Galactic) abundances of high-mass and both a thermal and a non-thermal component. The former is related 

low-mass X-ray binaries (HMXBs, LMXBs). The abundance to the SN-powered outgoing galactic wind and is mostly relevant at 

of SNRs (both Type II and la explosions) was also consistently energies < 1 keV. The latter is due to Compton scattering of the SN- 

estimated. From the literature PR02 derived typical spectra for accelerated, radio-emitting relativistic electrons off the FIR and CMB 

radiation fields, and it arguably dominates the spectrum of SFGs at 

Send offprint requests to: M.P.; e-mail: persic@ts. astro, it energies of > 30 keV (see Persic & Rephaeli 2002, 2003). 
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Persic et al.: 2-10 keV HMXB luminosity as galactic SFR indicator 




Fig. 1. The integrated "stellar" emission in the template X- 
ray spectrum of a star-forming galaxy proposed by Persic & 
Rephaeli (2002). Two different cases are shown: "standard" 
{left) and "top-heavy" (right). Shown for the standard case 
in ascending order at 7 keV: SNRs, faint LMXBs, HMXBs, 
and bright LMXBs. On the right, in ascending order at 3 keV, 
only SNRs and HMXBs are shown. The assumed 0.5-50 keV 
luminosities and number abundances are: logLx = 37.7 for 
HMXBs (50 objects) and high-luminosity LMXBs (70 ob- 
jects), 37.0 for SNRs (20 objects), and 36.7 for low-luminosity 
LMXBs (130 objects). The spectral components are normal- 
ized in energy flux in the 0.5-50 keV band. The spectrum is 
absorbed through a HI column density N^i — 10^^ cm^^. 



The PR02 approach has the flexibihty of handhng also the 
extreme cases of 'no ongoing SF' and 'very high ongoing SF' 
by switching off the HMXB and SNR components and, re- 
spectively, the LMXB component (see ¥ig.\-right). By letting 
the amplitudes of the various spectral components vary (while 
keeping their profiles fixed) in the spectral fit, one can use the 
PR02 procedure to determine the SF state of a given galaxy. 
Thus, this treatment is generally valid and not limited to any 
particular regime of SF activity. 

Based on their survey of galactic X-ray emission mech- 
anisms, PR02 concluded that in the 2-10 keV energy range 
X-ray binaries of both types, HMXBs and LMXBs, are the 
most prominent components with the required spectral shapes 
(see ¥\g.\-left). These have either a power-law (PL) form for 
HMXBs, or cut-off PL for LMXBs, with observed ranges of 
spectral parameters (HMXBs: F ~ 1.0-1.4, see White et 
al. 1983; LMXBs: F ~ - 1.4, ~ 5 - 10 keV, see 
Christian & Swank 1997) that provide good fits to the spatially 
unresolved fieppo5AZ spectra of the most extensively observed 
nearby starburst galaxies (SBGs), M 82 and NGC 253 (see 
PR02), and to the XMM-Newton spectra of a number of distant 
Ultra-Luminous Infra-Red Galaxies (ULIRGs; Franceschini et 
al. 2003). 

From a broader perspective, the template spectrum of PR02 
provides a framework that may prove especially useful for in- 
terpreting low spatial resolution data on SFGs, either distant 
(from Chandra and XMM: e.g., Hornschemeier et al. 2001, 
2003; Alexander et al. 2002; Bauer et al. 2002; Franceschini 
et al. 2003) or nearby (from A5CA and BeppoSAX: e.g., Cappi 
et al. 1999; Dahlem et al. 1998; Delia Ceca 1996, 1997, 1999; 
Moran et al. 1999). In particular, this template spectrum would 
allow the general possibility of measuring the ongoing star for- 
mation rate (SFR) in galaxies from their X-ray spectra, or per- 
haps - for some galaxies - directly from their X-ray luminosi- 



ties. The basic notion is that ongoing SFR in a galaxy can be 
measured based on stellar X-ray sources which are both suffi- 
ciently bright for their collective emission to be unambiguously 
identified, and sufficiently short-Uved so that they trace the 'in- 
stantaneous' SFR. Of the three main types of stellar galactic 
X-ray sources, 

(i) SNRs (which are X-ray bright over timescales ~ 10^ 
yr) are the evolutionary outcome of massive (8^ M/Mq^ 40) 
stars that explode on timescales (5^ r^v/Myr^ 50, see 
Maeder & Meynet 1989) short compared with a typical SB du- 
ration (rsB~ 100 Myr). Hence SNRs do trace the istantaneous 
SFR; however, their relative emission in integrated SFG spectra 
is quite modest and hard to identify; 

(ii) LMXBs (tx ~ 10^ yr) do contribute significantly to the 
X-ray emission but, due to the long delay between their forma- 
tion and the onset of their X-ray emission (the donor star has 

Mq), they do not trace the current SFR; finally, 
HMXBs (tx ^ 10^ yr) provide a suitable combination of 
short delay between binary formation and onset of X-ray emis- 
sion (the donor star has 8 Mq) and significant - some- 
times dominant - relative X-ray emission. 
Consequently, the measurement of ongoing SFR in galaxies 
hinges on our ability to separate out the HMXB contribution 
to the 2-10 keV luminosity, ^2-10- 

In this paper we extend the work of PR02 by examining 
ways in which L2-10 can be used as an astrophysicaUy moti- 
vated SFR estimator for SFGs. The possibility of using L2-10 
as an integral measure of the SFR has already been suggested 
(e.g.: Bauer et al. 2002; Grimm et al. 2003; Ranalli et al. 2003; 
Gilfanov et al. 2004). In this paper we propose that by using 
the HMXB portion of the hard X-ray luminosity, the SFR- 
L™^^ relation is universal and extends from normal galax- 
ies (low SFR) to very actively starbursting galaxies (very high 
SFR). Given the ways in which L2-10 can be contaminated, 

1. e. systematically from LMXBs at low-SFR regimes (PR02), 
and occasionally but quite frequently from AGNs at high-SFR 
regimes (Franceschini et al. 2003), our results will emphasize 
the potential and the limitations of using - over ~5 decades in 
X-ray luminosity - the 2-10 keV luminosity as an independent 
gauge of SFR in SFGs. 

The paper is organized as follows. In section 2 we discuss 
galactic SFR indicators linked to the presence of short-lived 
massive (^ SM©) main-sequence stars. Section 3 describes the 
sample of SFGs, which comprises objects with SBs of various 
strengths. In section 4 we discuss the effectiveness of the 2- 
10 keV luminosity of HMXBs as a SFR indicator. Section 5 
summarizes our main results. The conclusions are in section 6. 

2. Tracers of ongoing SFR 

As mentioned above, the indicators of current SFR are stel- 
lar objects or systems whose activity hfetimes are very short 
compared to characteristic SF timescales. Based on their spec- 
tra, luminosity, and number, HMXBs are most appropriate SFR 
gaugers in the 2-10 keV band. In these stellar systems the 
donor is a massive (M^ 8 Mq) OB star. Now, as a popula- 
tion OB stars are not important contributors to the X-ray output 
of galaxies (e.g., PR02), but their copious UV emission heats 
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Table 1. Data I: The sample of ultra-luminous IR galaxies (ULIRGs)^''). 



Object 




r 

/0.5-2 


r 

72-10 


Instr. 


P 

J60 


/lOO 


JNotes^ ' 




LMpcj 


LIO erg s J 


LIO erg s J 




UyJ 


UyJ 




IRAS 05189-2524 


170^1 




360.0'^1 


SAX 


13.25'**1 


11.84'*^' 


agn'^'I 


IRAS 12112-1-0305 


291 [3] 


L53P1 


1.62P1 


XIvIIvl 


8.18'^! 


9.46'=^1 


SB 


Mkn231 


169^1 


29.06P'"] 


84.5412.111 


XMM 


32.0'^1 


30.3'^1 


SB + AGN'21 


Mkn 273 


15l'^l 


17.00'*' 


70.0W 


ASCA 


2 1.7 [3] 


2i.4[3] 


SB + AGN'^1 


IRAS 14348-1447 


330l^I 


3.08^1 


1.91P1 


XMM 


6.82'=^! 


7.3l'^l 


SB 


IRAS 15250-3609 


2i3[31 


2.31^1 


2.31 '21 


XMM 


7.10'^1 


5.93'^' 


SB 


Arp 220 


73' ' 


8.0'''' 


18.0'"' 


SAX 


104.0'^1 


m.o'"'' 


SB 


NGC 6240 


97 [3] 


64.0[''l 


190.o[5.9l 


ASCA, SAX 


22.94'3l 


26.49 '^1 


agn"" 


IRAS 17208-0014 


170P1 


7.56P1 


4.54'21 


XMM 


9.53'31 


11.05'31 


SB 


IRAS 19254-7245 


246^1 


10.9il2.io] 


17.58'2'i°l 


XMM 


5.5'31 


5.8'31 


SB + AGN'21 
















(The Superantennae) 


IRAS 20100-4156 


517^1 


1.68^1 


2.24'21 


XMM 


5.2'31 


5.2'^1 


SB 


IRAS 20551-4250 


171(3] 


52.93 '21 


99.56'21 


XMM 


12.8'^1 


lO.O'^l 


SB + AGN'21 


IRAS 22491-1808 


309[3] 


0.61^1 


0.65'^1 


XIVIM 


5.54'31 


4.64'31 


SB 


IRAS 23060+0505 


692'^' 


270.00'^' 


250.0'^1 


ASCA 


1,2[31 


0.8''^' 


agn'^ 


IRAS 23128-5919 


17813] 


12.04'21 


21.66'21 


XMM 


10.8'^' 


ll.O'^^l 


SB -1- AGN'21 



References: [1] Brandt et al. 1997; [2] Franceschini et al. 2003; [3] Genzel et al. 1998; [4] Iwasawa 1999; [5] Iwasawa & Comastri 1998; 
[6] Iwasawa et al. 2001; [7] Severgnini et al. 2001; [8] Sanders et al. 2003; [9] Vignati et al. 1999; [10] Braito et al. 2003; [11] Braito et al. 
2004. 

'''^ Distances are taken from Genzel et al. (1998), or are computed consistently assuming Ho = 75 km s~^ Mpc~^, qo = 0.5 otherwise. 
Component(s) dominating the 2-10 keV emission; other name(s). 



up dusty molecular clouds, with the heat emitted in the FIR. 
Thus, the 2-10 keV emission of HMXBs and the FIR emis- 
sion of warm dust clouds are both related to the presence of 
short-lived OB stars (either isolated or in interacting pairs), 
with the emission ceasing soon after their death. In principle, 
then, the clouds' FIR luminosity and the 2-10 keV luminosity 
of HMXBs can be used independently to measure the instanta- 
neous SFR. The crucial step is the separation of the cloud FIR 
luminosity and the HMXB X-ray luminosity from the corre- 
sponding total luminosities: if this can be done properly, the 
FTR-derived and the X-derived estimates of the SFR should 
agree. 

Massive OB stars found in HMXBs usually end up as 
supernovae. Shock waves triggered by SN explosions ef- 
ficiently accelerate electrons to relativistic energies (e.g., 
Blandford & Ostriker 1980). Electron energy losses are radia- 
tive synchrotron-Compton at energies higher than ~ 300 MeV, 
with the relative significance of each of these processes de- 
pending on the value of the mean magnetic field and the scat- 
tering radiation field (either FIR or CMB; see, e.g., Rephaeh et 
al. 1991). Typical radiative energy loss times are t\oss ~ 10^ yr 
(e.g., Rephaeli 1979), i.e., comparable to or longer than a typ- 
ical SB duration. Therefore, galactic radio (synchrotron) emis- 
sion provides a measure of the SFR averaged over similarly 
long timescale, and has little to do with the ongoing SFR ex- 
cept in galaxies with constant SFR (see section 3). In the fol- 
lowing we shall deal with 2-10 keV luminosity of HMXBs and 
the FIR luminosity of dust clouds, both of which are immediate 
outcomes (and hence tracers) of the current SF activity. 



2.1. Infrared 

In normal disk galaxies, the relationship between the FIR lu- 
minosity and the SFR is complex because stars with a vari- 
ety of ages can contribute to dust heating, and only a frac- 
tion of the bolometric luminosity of the young stellar popu- 
lation is absorbed by dust (e.g., Lonsdale Persson & Helou 
1987; Walterbos & Greenawalt 1996). In SFGs studied here, 
the physical coupUng between the SFR and the IR luminosity 
is much more direct. Young stars dominate the radiation field 
that heats the dust, and the dust optical depths are so large that 
almost all of the bolometric luminosity of the SB is reradiated 
in the infrared. This makes it possible to derive a reasonable 
quantitative measure of the SFR from the FIR luminosity. 

A popular calibration of the SFR-Lfir conversion 
(Kennicutt 1998) is based on the SB synthesis models of 
Leitherer & Heckman (1995), which trace the temporal evolu- 
tion of the bolometric luminosity for a given SFR, metal abun- 
dance, and IMF. Kennicutt (1998) computed the SFR cahbra- 
tion using their "continuous SF" models, in which the SFR is 
presumed to remain constant over the lifetime of the burst. The 
models show that the Lboi/SFR ratio evolves relatively slowly 
between ages of 10 and 100 Myr, the relevant range for most 
of SBs (e.g., Bernlohr 1993; Engelbracht 1997). Adopting the 
mean luminosity for the 10-100 Myr continuous bursts, solar 
abundances, and the Salpeter (1955) stellar initial mass func- 
tion (IMF) defined as dN/dM oc M^'^-'^^ in the mass range 
0.1-100 Mq, and assuming the bolometric luminosity to be re- 
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Table 2. Data II: The sample of local normal and starburst galaxies (SBG sample)^°\ 



Object^'') 




Ad) 
JO. 5-2 


/2-10 


Instr. 




/so 


/lOO 


Notes'-'') 




[Mpc] 


[10-^^ erg/s] 


[10-^2 erg/s] 






[Jy] 


[Jy] 






1 o[261 
I. J 


1 q[2.20] 

1.0 


U.Oo 


POQ AT AQr'A 


/ .Oj 


77 nn[^] 


1 lA Oq[8] 

1 /4.uy 




isinr' no ^"5*^^ 


J.\J 


9 <[1.20] 


< n[i'201 


C AV 


7 no 


QA7 ft 1 [8] 

yo/.oi 


1 9S28 1 ^[8] 
IZoo.lJ 




iNO*^ UOZo 


Q 7(261 




fi 9/Iq[15] 

u.z4y 


A cr^ A 


7A 


91 Szl[^] 


szl 4*; [8] 

J4.4J ' 






Q ^[261 
y.O 


n s^[201 

U.<5j 


1 q[20] 

1.7 


A Qr^ A 


y.j 1 


AA /1a[^] 
00.40 


1 79 9^[8] 
1 /Z.Zj 






1.0 




fi 99 [3] 


A Qf^ A 


/19 


J4.J0 


cc 7q[8] 

jj.zy 




loUo 


lU.O 


U.UJ 


u. /o 


Aol_-/\ 




10s SS[^] 
iUJ.JJ 


141 7A[^] 






17 9(26] 
1 / .z 


U.OZ' 


1 Ills] 


A A 


10 SR 


1J.A AQ[8] 


1 QA OS [8] 
iy4.uj 




iN LiL. ZZ / 


Ofi o[26] 
JD.O 


U.Zl 


U.44 


A Qr^ A 
AoCA 


1 i. /J 


1 /I 9q[^] 

14. zy 


9C Q7[8] 

Zo.y/ 




NGC 2403^ 




1 A [20] 
1.0 


n Q^[20] 


A cr^ A 


0.4 J 


41.4/'^ ^ 


QQ 1 "3 [8] 

yy.ij"^ ^ 




iNOL^ Z /oZ 


07 ,[26] 


L.J 




POQ AT 


1 9 01 


17[8] 

y. 1 / 


1 'X 7 A [8] 
1 J. /o 






O.J 


7Q[20] 
u. 


U.DoD 


A <;p A 


Oil 


AO Szl[^] 


1 ^0 A^[^] 






r 9 [26] 


^ 8[21] 


IS ^[21] 




S SQ 
J.Jo 


1 AQC\ Al\-^] 
14oU.4Z ' 


1 ^7^ AO [8] 

1 J / j.oy 


— A;f 89 

— IVl oZ 


'wnr^ "5 mo 

JNvjL. DU/y 


on /i[26] 

ZU.'H-'^ ^ 


4 ';s[2]+ 

^.JO^ ^ 


n 7fi[2] 
yj. /o^ ' 


l?r*C AT A CP A 


1 y1 1 


'^O A7[8] 


1 C\A AO [8] 

iu4.oy 




INLJ*^ 3ZJD 


07 ^[26] 




U.JoO 


A A 


0.04 


1 09 A^[8] 


I 1 A 1 [8] 






10 7[26] 
lo./ 




91 [27] 


A <vP A 




^A sa[^] 

J4.J0 


AC 1q[8] 

46. ly 




NGC 3367^^ 


43.6^81 


0.18P°1 


0,16[^o] 


ASCA 


11.92 


6.44[«] 


13.48'*' 






14 1 [26] 
14.1 


U.44 


n ^n[20] 


A CP A 


y.O J 


^9 ^^[8] 

OZ.JJ 


7A Qo[8] 

/o.yu 


— 'wn 1 n8 

— iVl lUo 


L\KJ\^ JOZ.0 


7 7(26] 


4 1fil2] + 
H-.IO 


n Qs[2] 


DOC AT A CP A 


7. Ji 


'^A 80[8] 
J4.oU 


1 OS 7 A [8] 
lUJ. /D 




/\rp Z77 


41 filial 


U.J / 


1 /^o [20,23] 


A CP A 


1 1 RS 
i l.oJ 


ilJ.UJ 


1 1 1 d.9 [8] 


— INUl^ j\yy\j + Iv- D74 
















Anxr ot > in v*»\7[32] 
ALiiN at ^ lu Ke V 






79 [20,23] 


CO [20,23] 


A CPA 


1 n A9 

lU.OZ 


4J.10 


S7 00 [8] 

o/.uy 


= The Antennae 






U.oZO 


n zis9[4] 

U.'+oZ 


A CP A 


y.y4 


^A 0[11] 


ij.yj 








Ar\ n[2]+ 


1 1S[2] 
i.ij 


PMC AT A CP A 


S Al 
o.Dl 


Ao[8] 
oJ.4U 


1 AO OS [8] 
lOu.uo 






Ifi 8[26] 


n HA [20] 


r» no [20] 


A CP A 
AoCA 


lU. / J 


1 "5 "50 [8] 


'xn 77 [8] 
J 1 .1 r ^ 




NGC 4666 


14 l[26] 


0.16'"1 




SAX 


10.68 


37.1 iM 


85.95'*' 




NGC 4945° 


5.2P«1 


1.3[i°l 




SAX 


7.43 


625.46'*' 


1329.70'*' 


AGN 

at>10keV'l°'31.32,22,24] 


NGC 5236*"" 


47[26] 


3.5118] 


47[i8] 


ASCA 


7.98 


265.84'*' 


524.09'*' 


= M83 


NGC 5253° 




o'32'"l 




ROSAT 


10.47 


29.84'*' 


30.08'*' 




NGC 5457^^ 


5.4(26] 


0.54[^"1 


0,68['°] 


ASCA 


8.21 


88.04'*' 


252.84'*' 


= M101 


NGC 6946^ 


5.5[^«1 


3_0[20] 


1_2[20] 


ASCA 


7.78 


129.78'*' 


290.69'*' 




NGC 7469° 


65.2'"1 






ASCA 


12.64 


27.33'*' 


35.16'*' 


^Qj^[28,29,3a] 


NGC 7552° 


19.5I26] 


i.'oi^^i 




Einstein 


11.13 


77.37'*' 


102.92'*' 




NGC 7679 


71.0[^1 


3_3[6I 


6.0[^] 


SAX 


12.89 


7.40'*' 


10.71'*' 


agn'®' 


IC 342^^ 


39P6I 


L8P°1 


11.0[20] 


ASCA 


6.04 


180.80'*' 


391.66'*' 





References: [1] Cappi et al. 1999; [2] Dahlem et al. 1998; [3] Delia Ceca et al. 1996; [4] Delia Ceca et al. 1997; [5] Delia Ceca et al. 1999; 
[6] Delia Ceca et al. 2001; [7] de Naray et al. 2000; [8] Sanders et al. 2003; [9] Guainazzi et al. 1994; [10] Guainazzi et al. 2000; [11] Hunter et 
al. 1986; [12] Iwasawa et al. 1993; [13] Maccacaro & Perola 1981; [14] Martin & Kennicutt 1995; [15] Mizuno et al. 1998; [16] Moran et al. 
1999; [17] NED; [18] Okada et al. 1990; [19] Persic et al. 2003; [20] Ranalli et al. 2003; [21] Rephaeli & Gruber 2002; [22] Done et al. 2003; 
[23] Sansom et al. 1996; [24] Schurch et al. 2002; [25] Schulz et al. 1998; [26] Tully 1988; [27] Zezas et al. 1998; [28] Perez-Olea & Colina 
1996; [29] Nandra et al. 2000; [30] Blustin et al. 2003; [31] Madejski et al. 2000; [32] Delia Ceca et al. 2002. 

The superscripts G, R indicate whether an object is included in the Genzel et al. (1998) or Ranalli et al. (2003) samples, respectively. 
'■°) Distances are taken from Tully (1988) if cz < 3000 km s"'^, or are computed consistently assuming Ho = 75 km s"'^ Mpc~^ otherwise. 

Ruxes marked with a cross refer to the 0. 1-2.0 keV band. 
^"^ Blue apparent magnitudes, corrected to face-on and for Galactic absorption, from RC3 (de Vaucouleurs et al. 1991). 

Other name(s); spectrally dominating component(s) at energies ^ 10 keV. 

This lies within the range of published calibrations (1 — 3 x 
10"^"i\f : yr"^ (e.g.: Lehnert & Heckman 1996; IVleurer et al. 



10 '-^ (2.58 /eo + /loo) erg s ^ cm \ where /so and /loo are ex- 
pressed in Jy. 



radiated by dust, yields (Kennicutt 1998): 

^ The FIR flux is defined (Helou et al. 1985) as a combination of 
the IRAS 60/im and 100/im fluxes according to /fir = 1.26 x 



Table 3. Data III: The sample of Hubble Deep Field North galaxies (HDFN sample) 
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Source 


z 


FlAGHz 


I'1.4GHz 


SFrC-) 


/2-10 


^2-10 


Instr. 






[MJy] 


[ergs-^Hz-^] 


[Mq yr-i] 


[erg s~^cm~^] 


[ergs"^] 




134 


0.456 


210 


1.03 X 10=*° 


25.8 


2.8 X lO"^*^ 


1.38 X 10"*^ 


Chandra 


136 


1.219 


180 


8.68 X 10^^ 


217.1 


1.9 X 10"^® 


9.17 X 10"*^ 


Chandra 


148 


0.078 


96 


1.16 X 10^* 


0.3 


< 4.1 X 10"^^ 


< 4.95 X 10^* 


Chandra 


188 


0.410 


83 


3.23 X 10^^ 


8.1 


5.8 X 10"" 


2.26 X 10*° 


Chandra 


194 


1.275 


60 


3.24 X 10^° 


80.9 


2.0 X 10"^^ 


1.08 X 10*^ 


Chandra 


246 


0.423 


36 


1.50 X 10^^ 


3.8 


7.5 X 10"" 


3.13 X lO"*" 


Chandra 


278 


0.232 


160 


1.84 X 10^® 


4.6 


1.6 X 10"^® 


1.84 X 10"° 


Chandra 



All data are taken, or derived, from Ranalli et al. (2003). Distant dependent quantities have been computed assuming Ho = 75 km s ^ and 
qo = 0.5. 

The SFR has been computed using Condon's (1992) formula SFR(> 5 M©) = Li.4ghz/(4 x 10^* erg s"^Hz-^) M0yr-\ which yields 
the SFR for stars in the mass range 5 < M/Mq < 100 (assuming a Salpeter-like stellar IMF with dN/dM oc M~^'®, and a non-thermal 
radio spectral index a = 0.8). 



1997; Devereux & Young 1991; see also Inoue et al. 2000 for 
an analytical derivation, and Silva et al. 1998). 



Kennicutt's (1998) relation in eq.(l) is valid if all the FIR 
emission is due to warm molecular clouds heated up by newly 
bom massive stars: i.e., eq.(l) properly apphes to star-forming 
regions. If in addition to the warm component there is a "cirrus" 
component, i.e. emission from (colder) interstellar medium 
(ISM) heated by the general galactic UV radiation field which 
is powered also by relatively old stars, then a systematic bias 
will affect any SFR estimate based on eq.(l). Any cirrus com- 
ponent should then be filtered out from the observed FIR lumi- 
nosity, in order to determine the warm emission arising from SF 
activity. This correction is especially important for objects that 
are not SB-dominated, such as local normal galaxies and SBGs 
(where the SB affects only the central disk region), for which 
the cirrus component is expected to be relatively important. In 
order to do so we follow David et al. (1992) in using the results 
of Devereux & Bales (1989): Assuming that the O stars respon- 
sible for FIR emission are (ultimately) hnked to particle accel- 
eration (via SN shocks) and (in turn) to nonthermal radio emis- 
sion, it would be predicted that the cirrus component should 
be directly proportional to Lb (with the latter a measure of the 
galactic stellar content), i.e. L|fj^ = Lfir — xL-q. Devereux 
& Bales (1989) find that if x = 0.14, then for their sample of 
237 galaxies the FIR-radio (i.e., L|fj^-I/i.4GHz) correlation is 
optimized and becomes linear. This implies that, on average, 
14% of the blue luminosity of a normal or starbursting galaxy 
is reradiated in the FIR band as cirrus emission: for this class of 
galaxies (see Table 2), the FIR luminosity to be used in eq.(l) 
is then ilfp (which we here still refer to as LpiR for simplic- 
ity). Anticipating results obtained later in this paper, we can 
state that allowing for cirrus emission, while theoretically re- 
quired and practically favorable (doing so only strengthens our 
results), amounts only to a minor correction. 



2.2. X-Rays 

Being able to ascertain whether the 2-10 keV band is domi- 
nated by HMXBs or LMXBs may yield substantial insight on 
the SF history of a galaxy. In a SB environment the formation 
of low-mass stars could be suppressed if SN blast waves of 
more rapidly forming massive stars disrupt the slowly forming 
less massive stars before these complete (or even reach) their 
Hayashi tracks. The ensuing stellar IMF would be top-heavy 
(see, e.g., Doane & Mathews 1993 and Rieke et al. 1993 for 
M 82), in which case no LMXBs would form, and the resulting 
binary population of the SB consists only of HMXBs. But even 
if the IMF were not top-heavy, one isolated SB episode would 
not trigger an increase of the LMXB population. In fact, the 
time required for the ^ 1 Mq optical companion in a LMXB 
system to evolve out of the main sequence and come into 
Roche-lobe contact (and hence start the X-ray bright phase) 
substantially exceeds a typical SB lifetime (tsb~ 10® yr). So, 
whatever is the mass range of the relevant stellar IMF, during 
one isolated SB episode there is time only for HMXBs (i.e., by 
accretion-powered X-ray pulsars) to form (in larger numbers 
for a top-heavy IMF). The resulting 2-10 keV synthetic spec- 
trum would then be dominated by HMXBs, and hence the SB 
spectral slope would be characterized by a flat photon index 
r ~ 1.2 (seeFig.l-n^fo).^ 

In the case of recurrent SF bursts, LMXBs could be impor- 
tant contributors in the 2-10 keV band. This may happen in, 
e.g., galaxies that are members of pairs with highly eccentric 
orbits, or that are found in crowded environments (e.g., com- 
pact groups: see Hickson et al. 1989; see also Delia Ceca et al. 
1997). In both cases, tidal interactions would be recurrent, so 

^ Spectral components that are flat in the 2-10 keV band can come 
from AGNs. Both a direct nuclear emission (with F ~ 2) that is 
heavily absorbed, or an emission that is reflected by a warm medium, 
are observed as relatively flat (T ~ 1.2) spectra (e.g., NCC 6240, 
NGC 4945, Arp 299: Vignati et al. 1999, Guainazzi et al. 2000, Delia 
Ceca et al. 2002). In general there will be no degeneracy between the 
AGN and HMXB interpretations of the flat spectral components, as 
the inferred luminosity should enable a clear identification. 
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bright LMXB could be present as leftovers from earlier inter- 
actions. If one such galaxy is observed in a post-SB (or qui- 
escent) phase, the 2-10 keV emission would practically come 
only from LMXBs. 

Therefore, determining the type of X-ray binaries whose 
emission dominates the 2-10 keV luminosity of SFGs could 
provide a clue both to the understanding of SF history of those 
galaxies (e.g.. Holt et a. 2003), and to using L2-10 as an indi- 
cator of the ongoing SFR. As an example of the latter, let us 
consider three representative cases of SFGs. 

(i) In a galaxy where the SF activity has been constant over 
several 10^ yr, the relative frequencies of occurrence of stel- 
lar X-ray sources (SNRs, HMXBs, LMXBs) are expected to 
agree with those predicted for the Galaxy (which has a SFR of 
^ 2 Mq yr^^) by the Iben et al. (1995a,b) models, upon which 
the PR02 template X-ray spectrum is based. Comparing the de- 
duced HMXB content of a given galaxy (from X-ray spectral 
decomposition) with that of the Galaxy, its current SFR can be 
estimated. 

(ii) For SB-dominated galaxies the above considerations sug- 
gest that ^2-10 is an indicator of SFR. 

( Hi) In the intermediate case of a galaxy undergoing a minor 
SB (which is usually located in the central region), L2-10 is a 
superposition of the emissions from the (quiescent) bulge/disk 
and the ongoing SB. Only after these two components are sep- 
arated out, can the 'quiescent' and the 'bursting' SFRs be esti- 
mated (see Fig.2). 

3. The samples 

To test the idea that the SFR-Lx relation largely reflects the 
ability to use HMXBs to trace the current SFR, we have se- 
lected three representative samples of SFGs. 

The sample of ULIRGs (Table 1) - sources with bolomet- 
ric luminosity > emitted mostly in the IR (S-lOOO/xm) 

band - is essentially the list of IRAS-selected galaxies that rep- 
resent some of the most luminous galaxies in the local Universe 
observed with ISO by Genzel et al. (1998), whose FIR selec- 
tion makes it unbiased with respect to absorption. It is flux- 
limited and complete down to S'eo^im > 5.4 Jy, and high- 
quality IR/optical spectroscopic data are available for all the 
sources (Lutz et al. 1999; Veilleux et al. 1999). X-ray data in 
both the soft and hard bands are from ASCA, BeppoSAX, and 
XMM-Newton. According to the most extensive spectral sur- 
vey of ULIRGs to date (Franceschini et al. 2003), the ULIRG 
spectra are qualitatively similar to SBG spectra, plus possibly 
(in 40% cases) a heavily absorbed PL component that shows 
up at >5 keV (rest-frame). A similar incidence rate of AGNs 
in ULIRGs is deduced using mid-IR diagnostics (Genzel et al. 
1998). 

The sample of SBGs (Table 2), compiled from Genzel et 
al. (1998) and Ranalh et al. (2003), consists of local FIR- 
bright objects known to contain central starbursts of various 
strengths for which X-ray spectral information is available 
(from Einstein, ROSAT, ASCA, BeppoSAX, and RossiXTE). It 
should be stressed that this sample comprises a fair number of 
'normal' galaxies, i.e. galaxies whose SFR is comparable to 
that of the Milky Way. These SFGs are useful for checking the 



SFR-Lx relation at the low-SFR end of the galaxy distribution. 
A more appropriate name for this sample should then be 'lo- 
cal SFGs' (of which SBGs are, technically speaking, a subset); 
but for lack of a definite, clear-cut distinction between 'nor- 
mal' and 'starburst' galaxies, in this paper we'll stick to the de- 
nomination 'SBG sample' in order to emphasize that our main 
interest for these galaxies lies in their SF activity. Judging by 
currently available X-ray measurements, most of these galax- 
ies share a remarkable spectral homogeneity with their spectra 
usually described (most simply) as a soft thermal component at 
< 1 keV, plus a hard (sometimes cutoff) PL at > 2 keV (e.g., 
Dahlemetal. 1998). 

The sample of Hubble Deep Field North (HDFN) galax- 
ies (Table 3) was selected by Ranalli et al. (2003) from the 1 
Ms Chandra (Brandt et al. 2001) and radio (1.4 GHz: Richards 
2000 and Garrett et al. 2000) catalogues of the HDFN, which 
reaches limiting fluxes low enough for SFGs at z ~ 1.3 to be 
detected. As emphasized by Ranalli et al. (see also Grimm et 
al. 2003), these galaxies can be used to investigate the high-z 
behaviour of the SFR-X-ray luminosity relation, which is be- 
ing explored primarily in the local universe (z^ 0.1). Lacking 
FIR data, we estimate the SFR from the radio flux density at 
1.4 GHz using the relation SFR(> SMg) = ii.4GHz/(4 x 
lO^^erg s~iHz-i)MQyr-i (Condon 1992), wifli flie result- 
ing SFR are reported in Table 3. It should be emphasized that 
the Condon formula yields the formation rate of stars with 
5 < M/Mq < 100. To compare these SFR estimates with 
those we use for the SBG and ULIRG samples, we must ex- 
trapolate the radio SFR down to 0.1 Mq using the Salpeter 
(1955) stellar IMF upon which the Kennicutt (1998) formula 
is based. Doing so the radio and FIR estimates wiU be consis- 
tent ^ We flien obtain SFR(> 0.1 Mq) = 5.5 SFR(> 5 Mq). 
We will use these new values, SFR(> 0.1 Mq), as the refer- 
ence SFR values for our HDFN galaxies. It should be pointed 
out that, owing to the long Ufetime of the radio-emitting rela- 
tivistic electrons, ^^10^ yr, radio emission may not be able to 
trace shorter-timescale SBs (see section 2). However, it is re- 
aUstic to assume that HDFN galaxies, being in an earlier evo- 
lutionary phase than local galaxies, are dominated by strong 
continuous SF activity (e.g.: Madau et al. 1996; Thompson et 
al. 2001) and hence their SFR can be effectively measured by 
the non-thermal radio emission. 

4. HMXBs as gauges of the SFR 

We now demonstrate that the HMXB portion of the 2-10 keV 
luminosity, L™o^, can be used as a gauge of the SFR. To 
this end we compare the SFR computed from the FIR accord- 
ing to the Kennicutt (1988) formula [see eq.(l)] with the SFR 
estimated from the integrated luminosity of the bright-HMXB 
population. We assume a mean luminosity of 5 x 10^^ erg s~^ 
for bright HMXBs. For each sample galaxy, we estimate the 
number of HMXBs from the Lf^Q^. By comparison with the 
Galaxy, which hosts ~ 50 bright HMXBs (Iben et al. 1995a, 

In the derivation of the SFR-I/i.4,ghz conversion, Condon (1992) 
used a Salpeter-like IMF with index 2.50, instead of 2.35 as adopted by 
Kennicutt (1998). For our purposes here this difference is negligible. 
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M82: best 'Multi-kT Plasma + HMXB + LMXB + SNR' model 
ASCA/GIS+SIS + BeppoSAX/LECS+MECS+PDS 




chonnel energy (keV) 



M82: 'Multi-kT Plasma + LMXB + HMXB + SNR' madel 
ASCA + BeppoSAX data 




0.1 1 10 

channel energy (keV) 



Fig. 2. The "standard" model shown in Fig. 1 -left is fitted to the 
combined ASCA and BeppoSAX data of M 82 (ASCA: black, 
red, green, and dark-blue points; BeppoSAX: purple, cyan, and 
yellow points). The profiles of the spectral components are kept 
frozen in the fitting procedure (while their amplitudes are left 
free to vary) and {bottom) correspond to: photon index F = 
1.2 for HMXBs (yellow); photon index F = 1.4 and cutoff 
energy = 7.5 keV for LMXBs (purple); temperature kT — 2 
keV and chemical abundance Z — Zq for SNRs (cyan). At 
lower energies, three plasma components (having kT — 0.065, 
0.45, 0.75 keV; and Z = 0.1 Zq) are required (red, green, and 
dark blue). Shown are the model fitted to the data (top, energy 
spectrum), the residuals of the fit (middle), and the unfolded 
energy spectrum (bottom). 

and references therein), and has SFR ^ 2Mq yr^^ (e-g-, 
Matteucci 2002), we estimate the current SFR of the galaxy in 
question. This simple procedure is applied to both the ULIRG 
and SBG samples. 

• ULIRGs. If ^2-10 is used, the X-ray estimates are in overall 
agreement with the FIR estimates for the AGN-free ULIRGs 
(Fig. 4: filled dots), while they are clearly in excess for objects 
whose 2-10 keV emission is contaminated, or dominated, by 
a central AGN ^ (Fig. 4: empty dots and crosses, respectively). 

^ With both SFR estimates hnear in the respective luminosities, an 
equivalent statement is that, for a given Lfir, AGN-contaminated 
ULIRGs are X-ray overluminous compared with SB-dominated ones 
of the same Lpir; see also Ptak et al. 2003. 



IRAS 19254-7245 IRAS 20551 -4250 




^ . . . . . . .1 ^ P . . 1.1. II . . 

0.5 1 2 5 1 10 

channel energy (keV) channel energy (keV) 



Fig. 3. The energy spectra of the four AGN-contaminated 
ULIRGs of the Franceschini et al. (2003) sample. The XMM- 
Newton data for IRAS 19254-7245, IRAS 20551-4250 and 
IRAS 23128-5919 span the -0.4-10 keV energy range, while 
the composite XMM-Newton + BeppoSAX data for Mkn 23 1 
span the 0.2-60 keV energy range. The HMXB contributions 
have been estimated by perturbing the published best-fit spec- 
tral models with an additional F = 1.2 PL. Such additional 
components are reported in the four panels (clockwise from 
top left) as green, red, green, and red hnes. 



The overall agreement between FIR-derived and X-derived es- 
timates of the SFR suggests that the 2-10 keV luminosity of 
AGN-free ULIRGs is indeed (mainly) due to HMXB: this con- 
clusion is in accord with their flat observed 2-10 keV continua, 
which are consistent with the F ^ 1.2 spectra observed in 
Galactic HMXBs. 

For the objects in which the SB and AGN contributions 
to the 2-10 keV flux are comparable (see Table 1), we em- 
ployed the following crude procedure in order to estimate the 
HMXB contribution, a procedure that is admittedly quite ar- 
bitrary given the uncertain observational situation, but that 
nonetheless has some theoretical justification. For the 4 objects 
taken from Franceschini et al. (2003), we have added a F ~ 1.2 
PL component to the published spectral models ^, repeating the 

minimization. Results of the fitting, while not compelling 
because the quality of the data allows only simple model test- 
ing, do indicate a roughly ^ 10% HMXBs contribution to 
the 2-10 keV flux. Specifically, the additional PL fluxes are 
/2-10 = 2.06, 1.89, 1.44, and 9.25 x IQ-^^ erg cm'^ s-\ re- 
spectively (in the order listed in footnote 6), which in turn cor- 
respond to luminosities of L2-10 — 1-49 x lO''^, 6.59 x 10^°, 
5.44 X 10-^", and 3.15 x IG^^ erg s^^ (see Fig.3). For Mkn 273 
we have followed the argument of Xia et al. (2002) who, in 
their analysis of the Chandra spectrum, have identified a mod- 
erately absorbed PL component which they suggest must be 

" IRAS 19254-7245: Braito et al. 2003; IRAS 20551-4250, 
IRAS 23128-5919: Franceschini et al. 2003; and Mkn 231: Braito et 
al. 2004. 
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at least partially associated with HMXBs. (This PL is in ad- 
dition to a steeper and more absorbed PL that Xia et al. asso- 
ciate with the the direct light from the central AGN.) Therefore, 
lacking more detailed information, we assume the HMXB con- 
tribution in Mkn 273 to be half of the flux in the less-absorbed 
PL, i.e. /2-10 = 4.6 X 10~^^ erg cm^^ s~^ corresponding 
to -L2-10 = 1-27 X 10"^^ erg s~^. We emphasize that for the 
five AGN-contaminated ULIRGs of our sample the estimated 
HMXB contributions are very preliminary, given the substan- 
tial uncertainties. However, it is remarkable that the 'corrected' 
luminosities cluster so closely with the total 2-10 keV lumi- 
nosities of the AGN-free objects {Fig.5-right: empty dots). 

• SBGs. In comparing SBGs and ULIRGs, we should bear in 

mind a main difference between the two classes that is par- 
ticularly relevant to our discussion here. In ULIRGs the star 
forming activity is generally very intense and involves most 
of the stellar disks; in local SBGs (such as those included in 
Table 2), SF is considerably less intense and occurs mainly in 
the very central disk region. Thus, for non spatially resolved 
measurements, such as the ASCA and BeppoSAX data used in 
Figs.4 and 5, the 2- 10 keV emission is expected to come mostly 
from HMXBs in ULIRGs, and from a mix of HMXBs (in the 
central starbursting region) and LMXBs (in the underlying qui- 
escent disk) in normal and staburst galaxies. Therefore, if this 
effect is not corrected for, our X-ray-based estimate of the SFR, 
while fairly adequate in principle (albeit crude) in the case of 
ULIRGs, is systematically overestimated in the SBG sample by 
a factor 1//, with / the HMXB-to-total luminosity ratio. Thus, 
we need to estimate / for the SBG sample listed in Table 2. 

To this end we fit the template X-ray model of PR02 (see 
Fig.l-left), keeping the profiles of the spectral components 
frozen, and letting only the amplitudes free to vary, to available 
data for such galaxies. In principle this procedure should be ap- 
pUed to each and every galaxy in the sample. In practice, how- 
ever, for most of the objects the available data are such that it 
is difficult even to determine whether the harder spectral com- 
ponent is PL or thermal (e.g., Dahlem et al. 1998; Dahlem et 
al. 2000), let alone a more complicated multi-component (ther- 
mal and PL) model including also emisssion from HMXBs, 
SNRs, and LMXBs. The very nearby M82, which has been ob- 
served with ASCA in the 0.4-10 keV range (Ptak et al. 1997) 
and BeppoSAX in the 0.1-60 keV range (Cappi et al. 1999), is 
currently the best example of a relatively well-observed SBG. 
In Fig. 2 we show the result of fitting the PR02 template to the 
combined ASCA+BeppoSAX data (see the Appendix). The re- 
sults of these fits yield / ~ 0.2. 

In order to check to what extent this value of / may be 
taken as representative of SBGs in general, we have performed 
additional fits for a subset of our SBG objects (including M 83, 
NGC 253, NGC 2146, NGC 2903, NGC 3310, NGC 3256) be- 
ing aware, however, that the archival A5CA spectra available for 
this analysis are inferior to that of M82, so any result from such 
a comparison is necessarily preliminary and merely suggestive. 
(For the other objects of the SBG sample for which ASCA data 
exist, we felt that the quality of the data was not sufficient for 
even a rough estimation of /.) In addition to one or more sub- 
keV thermal plasma components, the fitting models we used 



did include a HMXB component plus either a plain/cutoff PL 
(NGC 2146, NGC 2903, NGC 3256, NGC 3310), or the ex- 
phcit SNRh-LMXB components (M 83, NGC 253). Wifliin the 
hmits imposed by the relatively poor statistics of these data, the 
results are consistent with / ~ 0.2 for the objects with high- 
est quahty data (M 83, NGC 253, NGC 2903), while suggest- 
ing a sUghtly higher value of / in the other cases (NGC 2146, 
NGC 3256). The exception seems to be NGC 3310, whose hard 
component is flat enough (F ~ 1.4: Zezas et al. 1998) to be 
accounted for by HMXBs alone, and hence / ^ 1. In con- 
clusion, our additional spectral analysis supports the deduction 
that SBGs have / ~ 0.2, while showing a high scatter in the 
actual value of this factor. 

The possibility that the emission of SBGs may be domi- 
nated by a few very luminous point sources (e.g., for M 82 
see Kaaret et al. 2001 and Strohmayer & Mushotzky 2003), 
and that analyses of spatially integrated SBG spectra may be 
biased, motivates further attempts to estimate /. Information 
about the relative emission from HMXBs can also be derived 
from studies on point source populations in nearby galaxies 
based on the high angular resolution Chandra data. Studies 
on individual local galaxies have yielded luminosity functions 
of X-ray point sources (XPLFs) down to Umiting luminosi- 
ties of -lO^^""^^ erg s-i (e.g., Kilgard et al. 2002; Colbert 
et al. 2003; Fabbiano et al. 2001; Zezas et al. 2002, 2003; 
Roberts et al. 2002a; Soria & Wu 2002; Bauer et al. 2001; 
Griffiths et al. 2000; Swartz et al. 2003; Kaaret 2002; Soria et 
al. 2003; Gao et al. 2003; Terashima & Wilson 2003; Kong 
2003). Such XPLFs can be described as N{> L) oc L"" 
with a ~ 0.5, 1, 1.5 for starburst, normal spiral, and ellipti- 
cal galaxies (and bulges), respectively (Kilgard et al. 2002 and 
references therein). The difference in XPLFs among galaxy 
types, which implies that higher-SFR galaxies have relatively 
more sources (presumably, HMXBs) at high luminosities, is 
presumably related to the age of the underlying X-ray binary 
population. A simple birth-death model of X-ray binaries pre- 
dicts a difference of 1 between the XPLF slope for an equilib- 
rium XP population and that for an impulsively-formed quietly- 
ageing XP population. However, XP counts do not give clues 
as to the nature of the XPs. Prestwich et al. (2002) have pro- 
posed to classify XPs in external galaxies on the basis of their 
X-ray colors (defined using Chandra's 0.65, 1.5, and 5 keV 
bands), finding highly significant differences in the X-ray col- 
ors of bulge and disk XPs, i.e. of LMXBs and of HMXBs and 
SNRs. Colbert et al. (2003) have analyzed 1441 XPs detected 
in 32 nearby galaxies, concluding that the XP luminosity is 
well correlated with the K-luminosity (i.e., the stellar mass) 
and the FIRh-UV luminosity (i.e., the SFR), suggesting that 
XPs are connected to both the old and young stellar popula- 
tions. Using their X-ray color diagrams, we find that ~0.2 of 
the XPs observed in spiral and interacting galaxies can be clas- 
sified as HMXBs. Correspondingly, a fraction ~0.25 of the 0.3- 
8 keV XP luminosity of these galaxies can be attributed to such 
F '--^ 1.2 sources. Because the XP luminosity typically accounts 
for ~ 80% - 90% of flie total luminosity (e.g., Roberts et al. 
2002a), we deduce that (putative) HMXBs account for ~20% 
of the integrated 0.3-8 keV luminosity. This value is in reason- 
able agreement with our spectral estimate of / ~ 0.2, mostly 
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based on M 82. ^ We therefore consider / ~ 0.2 to be the cur- 
rently appropriate, representative average value, for the SBG 
sample. 

An analogous correction to the FIR emission, aimed at re- 
moving any contribution that is unrelated to ongoing SF activ- 
ity, is also in order. As discussed in section 2.1, we must sub- 
tract the (SB-unrelated) cirrus component from the FIR emis- 
sion of SBGs ^. [The practical impact of this necessary correc- 
tion is minor (see Fig.4), and thus has no bearing on any of our 
results.] 

The SFR, as measured by /L2-10, turns out to be in 
good agreement with that inferred from Lfir (see Fig. 5, filled 
squares ^); the left panel includes the correction for the cirrus 
component in the FIR emission. To appreciate the effect of the 
correction, in the right panel we show the plot without the cir- 
rus correction. 

• Hubble Deep Field North Galaxies. In the SFR-i2-io plane 
the HDFN galaxies overlap (with a wider spread in SFR) with 
the SB-dominated ULIRGs (see Fig.5-/e^) i". 



The integrated XP luminosity of SFGs is dominated by ultra- 
luminous X-ray sources (ULXs), i.e. off-nucleus objects with L > 
10'^^ erg s^^ which is the Eddington luminosity for a 8 Mq black 
hole (BH) (which is the limiting BH mass obtainable via stellar evolu- 
tion). The nature of ULXs is still unclear. However, ULXs appear to be 
related mostly to high-SFR environments (Fabbiano et al. 2001; Lira 
et al. 2002; Roberts et al. 2002a; Zezas et al. 2002, 2003; Humphrey 
et al. 2003; Gao et al. 2003; Karet et al. 2001; Strickland et al. 2001): 
the identified optical counterparts of ULXs are indeed O stars, so sug- 
gesting a relation of ULXs with HMXBs (Liu et al. 2002; Roberts et 
al. 2002b). The spectra of ULXs, reminiscent of those of Galactic BH 
X-ray binary candidates in high state, are described by a model com- 
prising an accretion disk with inner-edge temperature of ~1.1 keV 
plus a hard (F ~ 1.2) PL similar to that observed in HMXBs (Zezas 
et al. 2002; Foschini et al. 2002; see also Fabbiano & White 2003). 
A scenario in which a population of HMXBs, with ^ 15 Mq donors 
transferring mass on their thermal timescales (and hence at super- 
Eddington rates) to ^ IOMq accretors, fits current observations on 
ULXs (King 2003). Our computation of / ~ 0.2, based on the XP 
colors corresponding to F ~ 1.2, does therefore incorporate the pos- 
sible presence of ULXs. 

Some objects in Table 2 are found at low Galactic latitudes 
(|6| < 15 degrees). For these the corrections for foreground Galactic 
absorption tend to be large and uncertain. An extreme case is IC 342 
(b — 10.58 deg), for which the B-band absorption is estimated to be as 
large as 3.360 mag (Burstein & Heiles 1982) or 2.407 mag (Schlegel 
et al. 1998). As the adopted statistical correction for cirrus emission 
turns out to be unphysical in this case, we will leave IC342 out from 
further analyses involving cirrus-corrected FIR luminosities. No re- 
sult in this paper will depend on whether this object is included in (or 
excluded from) the sample. [Note that IC 342 is uncharacteristically 
underluminous in FIR and soft X-rays compared with hard X-rays, 
suggesting the presence of a HMXB population typical of a SBG, but 
with reduced FIR and hot gas emissions (Bauer et al. 2003)]. 

' The SBG sample contains one object, Arp 299, which was re- 
cently shown to host an AGN (Delia Ceca et al. 2002). Remarkably, 
this object lies in the region populated by the AGN-dominated 
ULIRGs. 

Had we mistakenly used the plain SFR(> 5 Mq) values resulting 
from applying Condon's (1992) formula to the radio data - instead of 
the extrapolated SFR(> 0.1 Mq) values - the HDFN galaxies would 
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Fig.4. The SFR estimated from X-ray emission versus that 
estimated from FIR emission. For all objects it is assumed 
that the emission in both bands is related to starburst activity. 
Filled squares and filled dots represent, respectively, SBGs and 
SB-dominated ULIRGs; empty dots denote SB/AGN-powered 
ULIRGs, while crosses denote AGN-dominated objects; filled 
triangles denote HDFN galaxies (for which, due to the lack of 
FIR measurements, the SFRs have been estimated from radio 
fluxes). The dotted line, representing the SFRx = SFRfir re- 
lation, is only meant to guide the eye. The right panel shows 
the case where the FIR luminosities of normal and starbursting 
galaxies are not corrected for cirrus emission. 



Lacking spectral information, a simple and realistic as- 
sumption is that HDFN galaxies - likely 'juveniles' - are SF- 
dominated. Indeed Cohen (2003) finds that, based on the emit- 
ted luminosity in the 3727 A line of [O II], the average SFR is 
about an order of magnitude higher in HDFN galaxies than in 
local ones. Furthermore, the age of a galaxy at z ~ 1, ~6 Gyr, 
is not long enough for LMXBs - which have sub-solar mass 
donors - to form (see Maeder & Meynet 1989). These consider- 
ations suggest that the 2-10 keV emission of HDFN galaxies is 
due essentially to HMXBs, like in the SB-dominated ULIRGs 
discussed above. We then set L™^^ = 1/2- 10 • The X-ray es- 
timates of the SFR agree comfortably with those inferred from 
the radio data (see Fig.4). Therefore, the SFR-Lf^^^^^ relation 
for HDFN galaxies is consistent with that for SBGs and SB- 
dominated ULIRGs (see Fig.5-right). 

Analysis of the three samples leads us to the following con- 
clusions: 

/) In SB-dominated ULIRGs the 2-10 keV luminosity is a 
gauge of the SFR. This result lends further support to the in- 
terpretation of the ~ 2-10 keV emission of ULIRGs as being 
(mainly) the integrated emission from HMXBs. In those cases 
where an AGN clearly shows up (Franceschini et al. 2003) 
when this component is subtracted out, the residual 2-10 keV 
luminosities imply SFR values that agree with the FTR-based 
ones. 

i) In local star-forming (i.e., normal and starburst) galaxies only 
a relatively minor fraction (/ ^ 0.20) of the 2-10 keV luminos- 
ity is related to the ongoing SF activity. Correcting for this bias 
results in essentially the same SFR-L™j^^^ relation for both 

ULIRGs and SBGs. 



have lined up, and partially overlapped, with the SBGs in the SFR- 
L2-10 plane. 
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Fig. 5. The SFR versus 2-10 keV luminosity relations, using 
the total luminosity (/e//) and the HMXB luminosity (right). 
The SFR has been computed from the cirrus-subtracted FIR 
luminosities according to eq.(l) [except for the HDFN sam- 
ple (filled triangles) where radio data have been used; only X- 
ray detections have been used, see Table 3]. Symbols are as in 
Fig.4; the large empty square represents the Milky Way (SFR ~ 
2 Mq yr-\ see Matteucci 2002; L2-10 ~ 2 x 10^^ erg s"^ and 
Lf-w^ ~ 4x 10^^ ergs" \ estimated from Ibenetal. 1995a,b 
and consistent with the template spectrum in Fig.l). The solid 
line in the right panel shows the linear relation in eq.(2), 
whereas the dashed line in the left panel shows the Ranalli et 
al. (2003) relation SFR = L2-io/(5 x lO^^erg s-^)MQyi-'^. 

Hi) Distant (z ~ 1) SFGs have SFR values and 2-10 keV lumi- 
nosities similar to those of ULIRGs. 

Within the uncertainties of the data (most notably those in- 
volved in identifying the HMXB contribution in our sample 
of AGN-contaminated objects) the resulting relation is roughly 
linear, 



SFR(>O.1M0) = 



rHMXB 
-^2-10 

lO^^erg s~ 



Moyr" 



(2) 



with an estimated 20% statistical error. 

It should be remarked that the SFR-luminosity relation 
is not tight if no correction is applied (see Fig.5-Ze/f). Both 
SBGs and AGN-contaminated ULIRGs are overluminous for 
their SFR when compared with the SB-dominated ULIRGs and 
HDFN galaxies. Indeed, SBGs are contaminated by substantial 
LMXB emission, and in AGN-contaminated ULIRGs the AGN 
emission is typically brighter than the SB emission. After cor- 
recting for both types of excess, no hint of non-linearity appears 
in the SFR-luminosity relation (Fig.S-right). 

5. Discussion 

Our main suggestion - and the central theme of investigation 
- has been the notion that the 2-10 keV collective emission 
of HMXBs, i™i"o^, is a meaningful gauge of ongoing galac- 
tic SFR. As spelled out in Section 4, given the universaUty of 
the accretion process onto NS/BH in binaries, the linear depen- 
dence of LpiR on the SFR, the direct link between the FIR and 
HMXB emissions, a tight relation is predicted between SFR 
and L2^^^. Whether the latter relation is also linear depends 
largely on whether the stellar mass range and the shape of the 
IMF are the same in all galaxies, or vary systematically with 
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Fig. 6. The X-ray-FTR luminosity relation for soft [mostly 0.5- 
2 keV (left)] and hard X-rays [2-10 keV: total emission (cen- 
ter), and for just HMXB emission (right)]. In both panels the 
FIR luminosities of the SBG sample have been corrected for 
cirrus emission. Symbols are as in Fig.4. 

the SFR. In the former case we expect the relation to be lin- 
ear, whereas if the IMF becomes progressively more top-heavy 
(i.e., the lower mass cut-off becomes progressively higher) with 
increasing SFR, a non-linear relation is expected. 

As we have mentioned, the use of HMXBs as a galactic 
SFR estimator was already proposed by Grimm et al. (2003) 
(see also David et al. 1992). Within the general agreement be- 
tween their background picture and ours, the main result of 
Grimm et al. is, however, quite different from ours: their result- 
ing relation is non-linear (SFR oc L°'^) for small SFR and low 
values of L2-10 (i-e., SFR < 4.5 Mgyr"^ and ^2-10^ 2.6 x 
10^°erg s~^), and linear for higher values of the SFR and 
L2-10- Gilfanov et al. (2004) argue that such non-Unearity in 
the low-SFR limit may be caused by non-Gaussianity of the 
probability distribution of the integral distribution of discrete 
sources. We can offer no clear explanation for the discrepancy 
between our result and theirs. One important difference be- 
tween our approach and that of Grimm et al. (2003) is their 
a priori minimization of the LMXB contribution because - as 
they put it - 'owing to the absence of optical identifications of 
a donor star in the X-ray binaries detected by Chandra in other 
galaxies, ... there is no obvious way to discriminate the con- 
tribution of low-mass X-ray binaries' . Grimm et al. argue that, 
given the long evolution timescales of LMXBs, in a galaxy the 
LMXB population should be proportional to the galaxy stellar 
mass, whereas the HMXB population should be proportional to 
the SFR, so that the relative importance of LMXBs should be 
characterized by the ratio of the stellar mass to the SFR. Using 
dynamical estimates for the stellar masses of galaxies, and SFR 
derived from a variety of indicators (UV, Ha, FIR, radio, that 
give a wide range of values of the SFR - see their Table 3), 
Grimm et al. claimed that they selected only galaxies where 
HMXB emission is expected to exceed LMXB emission by a 
factor of ^ 3. For their sample galaxies, consequently, Grimm 
et al. supposed that the plain 2-10 keV luminosity is a measure 
of the collective HMXB emission in that band. In contrast, our 
procedure involves the identification of the HMXB emission 
from a given SFG by means of a spectral decomposition of the 
observed X-ray spectrum. 

The need in our approach to single out 0^ rather than 
L2-10 is further spurred by the mix of starbursts of widely 
varying strengths, from low/moderate (local SBGs) to extreme 
(ULIRGs). While the FIR emission is mostly related to cur- 



11 



§ 

5 




0.0001 0.001 0.01 0.1 
L(3-10 keV) / L(FIR) 



Fig. 7. The distribution of 60/i-to-lOO/i flux density ratios as 
a function of the SFR {left) and of the SB 2-10 keV to FIR 
luminosity ratio {right). Because /eo/Zioo is an approximate 
measure of the temperature, these plots show that: ( i) in spite 
of a relatively large scatter, a SFR-temperature correlation is 
present in the data, whereby galaxies with higher (lower) SFR 
tend to have higher (lower) dust temperature, T ^ 50 (25) °K; 
and (ii) the ratio of HMXB 2-10 keV luminosity to dust FIR 
emission remains roughly constant irrespective of the dust tem- 
perature. Symbols are as in Fig.4. 



rent SF in both cases (see section 2.1), the 2-10 emission 
is not in the former case but it is in the latter. Had we re- 
stricted to the local SBG sample, we would have found SFR 
— -^2-10/(5 X lO'^'^erg s~^)AfQyr^^ (see Fig.5-Ze/0, in agree- 
ment with Ranalli et al. (2003). In fact, the 2-10 keV emission 
from local SBGs has approximately a similar mix of SB-related 
(HMXB) emission and SB-unrelated (LMXB) emission, so that 
for this sample the primary correlation, SFR oc L^^q^, prop- 
agates into SFR oc L2-1Q. Some extra scatter plagues the latter 
relation because the value / ^ 0.2 appearing in the equality 
j^sMXB _ J ^2-10 most likely has some appreciable scat- 
ter within the sample. Only when the SB-dominated ULIRGs 
and HDFN galaxies are introduced into the plot and are used 
as calibrators, does the linearity of the SBG SFR-L2-10 re- 
lation break down (see Fig.5-left). The Unearity of the rela- 
tion, extended to the whole luminosity range, is re-established 
by using the SB-related luminosity i™^^ (see Fig.5-right). 
[Of course, the calibration of our relation is lower than that of 
Ranalli et al.'s relation by the factor / = 0.2.] 

One key assumption underlying our SFR-L2-10 relation 
is that the FIR luminosities - from which the SFR values are 
computed - are really due to the SB and are not contami- 
nated by other contributions in any important way. In order 
to check this assumption, we convert the observed /eo //loo 
ratios to dust temperatures using the table in Helou et al. 
(1988). The distribution of 60/x-to-lOO/i flux density ratios (see 
Fig.7 -left) implies a range of estimated dust temperatures 
25 T^S 50 °K (if the dust emissivity is proportional to 
i^~^), the lower bound typically referring to more quiescent 



Because galaxies have multiple emission components with widely 
varying parameters, the ratio /eo //loo cannot be expected to yield a 
precise temperature. The uncertainty affecting dust temperature esti- 
mates is largely intrinsic, and can be traced back to substantial uncer- 
tainties that are essentially of astrophysical origin (e.g., the frequency 
dependence of the dust emissivity, e.g., Helou et al. 1988). 



(i.e., 'normal') galaxies (e.g., NGC 628) and the upper bound 
to more actively star forming objects (e.g., the well-known SB- 
dominated galaxy Arp 220; Silva et al. 1998). 

This range of values seems realistic for the samples of SFGs 
considered in this paper, thus confirming the soundness of the 
assumption that our FIR luminosities are not significantly con- 
taminated by emission unrelated to ongoing SF activity, and 
verifying earlier suggestions that galaxies with higher SFRs 
have higher dust temperatures (e.g., David et al. 1992; but there 
may be exceptions to this trend, see FigJ-left). Notice also that 
the ratio of the two SB-related luminosities, L™^^/Lfir, re- 
mains approximately constant when moving from lower tem- 
peratures (normal and starburst galaxies) to higher tempera- 
tures (ULIRGs), log (La-io/ipm) ^ -4.3 ± 0.3 (see Fig.7- 
right). In particular, note that, when properly accounting for 
LMXB and AGN emission, SBGs and ULIRGs imply an X-to- 
FIR ratio consistent with that for the pure SBs arguably repre- 
sented by the AGN-free ULIRGs. 

In the simple scenario, in which the formation of OB stars 
and of accreting NSs/BHs during episodes of active SF activity 
only depends on local conditions, the same FIR-X-ray emis- 
sion relation is expected to hold for local (z^ 0.1) and distant 
{z ~ 1) star-forming galaxies. Indeed, the sample of 2; ~ 1 
HDFN galaxies we have used does follow the SFR-Lx rela- 
tion estabUshed in local galaxies. It should be emphasized that 
distant galaxies are in a substantially active SF phase: their 2- 
10 keV emission must then be HMXB-dominated, if no central 
AGN is contaminating the emission. If an independent estimate 
of the SFR is also available (e.g., from deep radio observa- 
tions), the location of the galaxy in the SFR-L2-10 plane can 
be used to evaluate any excess X-ray luminosity which, if de- 
tected, would presumably imply the presence of an AGN. For 
example, if according to eq.(2) a distant galaxy were very over- 
luminous for its SFR (e.g., if ^2-10 ~ 10^'^ erg s~^ for SFR 
~ 100 Mq yr~^), we might suspect a significant AGN contri- 
bution to the observed 2-10 keV emission (evolutionary con- 
siderations rule out significant LMXB emission). This method, 
based on joint X-ray and (say) radio observations, could be a 
very efficient tool to detect AGNs that are buried deeply (i.e., 
absorbed through iVn > 10^^ cm~^) in the nuclei of distant 
SFGs and are not detectable in other bands [e.g., the exemplary 
cases of NGC 6240 (Vignati et al. 1999), NGC 4945 (Guainazzi 
et al. 2000), Arp 299 (Delia Ceca et al. 2000)]. We therefore 
caution against using X2-10 as a gauge of the SFR. 

Finally, from Fig.6-(Ze/0 we see that the emission of SFGs 
at ~1 keV is also an indicator of the SFR, although the scat- 
ter of the SFR— Lx relation is higher in the soft band than in 
the hard band (see Fig. 6, where the equivalent ipiR-ix rela- 
tion is shown). Part of the scatter seen in Fig.6-{left) is due to 
observational errors - such as uncertainties in determining the 
absorption and metallicity of the soft thermal plasma, or occa- 
sional differences in the definitions of the soft spectral bands 
- but most of it is probably intrinsic. In fact, the soft compo- 
nent which is systematically observed in the spectra of SFGs 
is interpreted as sub-keV thermal emission originating from 
SN-powered outgoing galactic winds (e.g., Dahlem et al. 1998; 
Franceschini et al. 2003). Spatially resolved data (Strickland 
et al. 2000) and simulations (Strickland & Stevens 2000) have 
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shown that this emission occurs at the boundary between the 
hot, tenuous wind fluid proper and the cool, denser ISM where 
conditions are optimal for thermal emission of X-rays. This im- 
pUes that only a small fraction of the wind mass is involved in 
X-ray emission, and that local conditions (i.e., density, dumpi- 
ness, and chemical composition of the ISM) crucially deter- 
mine the level of this emission. The small mass fraction of ther- 
mally X-ray emitting galactic gas is subject to a large scatter 
among galaxies. These considerations suggest that most of the 
scatter observed in Fig.6-(/e//) is probably intrinsic. Indeed, the 
pure thermal emission does not generally correlate with Lfir 
(see Fig. 6c of Franceschini et al. 2003), implying that galactic- 
wind emission is not a SFR indicator (in spite of winds being 
an immediate SB outcome). Then, we suggest that the observed 
^0.5-2 — i'FiR correlation is not primary and mainly due to 
SB-powered galactic winds in the SFGs, but that it is probably 
induced by the portion of the hard spectral component of SFGs 
showing up in the 0.5-2 keV band. Notice how the Lx — ^fir 
relation improves when going from left to right in Fig.6, i.e. 
when the reference X-ray luminosity changes from 'total soft' 
through 'total hard' to 'HMXB hard'. 

6. Conclusion 

Current hmited ability to spectraUy decompose the 2-10 keV 
emission of SFGs is nonetheless indicative of the feasibility 
of using the main HMXB component as an estimator of the 
galactic SFR. The total 2-10 keV luminosity is not a proper 
SFR indicator since for low values of the SFR - in normal and 
moderately-starbursting galaxies - it is substantially affected 
by the emission of LMXBs, while for high values of the SFR 
- in ULIRGs - it usually includes significant emission from 
obscured AGNs. 

The use of an additional sample of distant galaxies from the 
Hubble Deep Field North, for whom there is no X-ray spectral 
information, but whose SFR can be estimated from deep radio 
data, led us to to conjecture that the SFR-Lx relation that we 
deduced for local galaxies might also hold for distant, ^ ~ 1, 
galaxies, if the 2-10 keV emission of HDFN galaxies is domi- 
nated by HMXBs. 

The SFR-i™^^ relation may be universal across galaxy 
types and redshifts. This may not be so surprising, if our 
understanding of the processes involved is basically correct. 
However, its linearity and caUbration, both in the local and in 
the distant Universe, are major issues that remain to be set- 
tled. Resolved high energy long-exposure XMM-Newton spec- 
tra will be needed to afford detailed modelling of the integrated 
emission, while high-spatial-resolution Chandra imaging will 
be needed to construct the luminosity functions of HMXBs 
Achieving mutually consistent results from the spectral de- 
composition and from the integration of the HMXB luminosity 
function would lead to a definite determination of L™™. 

Once the universality of the SFR-L™™ relation has been 
cleary established, it will be possible to use the relation as a tool 
to unveil the presence of deeply absorbed AGNs lurking in the 
centers of distant galaxies. If, for its (independently estimated) 
SFR, the 2-10 keV luminosity of a distant galaxy will be mea- 
sured to be in excess of that predicted based on the SFR-Lx 



relation, that galaxy would be suspected to harbor a strongly 
absorbed AGN. 

Finally, we stress again that in order to build a physically 
sound X-ray SFR indicator, we should use, in principle, a sam- 
ple of star-forming galaxies whose X-ray spectral properties 
are well known, and for which the SFR estimates are robust. 
Extracting the HMXB 2-10 keV emission is clearly one crucial 
step. Ideally, high-quality spectra should be available and ana- 
lyzed for each object in the sample. We unfortunately do not yet 
have such data. At moderate- and low-SFR regimes we relied 
on estimates gained from analyzing data for few objects (M 82, 
M 83, and NGC 253; plus educated guesses from NGC 2146, 
NGC 2903, NGC 3310, NGC 3256), which we then assumed 
would hold for the whole sample of normal and starburst galax- 
ies. In the high-SFR regime we do have suitable spectral anal- 
yses for the purely SB-powered ULIRGs, but we had to resort 
again to estimates in the case of AGN-contaminated ULIRGs. 
For the sample of distant HDFN star-forming galaxies we could 
only make (educated) assumptions. The situation is clearly 
far from ideal. Therefore our results, although promising are 
clearly preliminary. With the ongoing acquisition of spatially 
resolved Chandra data and spectrally resolved XMM-Newton 
data, it should be possible to observe galactic regions that are 
clearly recognized sites of ongoing SF, in particular the ac- 
tively star-forming central regions of starburst galaxies. These 
smaller starbursts, along with those ULIRGs which are not con- 
taminated by AGNs and can be considered as giant starbursts, 
could yield a well defined SFR-Lx relation. So far we have 
been able to suggest a continuity of SF properties moving from 
normal through starburst to ultra-luminous IR galaxies (quan- 
tified by L2-^^io^/-Z^FiR =i 5 X 10"^ over about 5 decades in 
X-ray luminosity). But the SFR-ix relation we have obtained 
is still affected by a fairly large scatter, which would currently 
hmit its applicabiUty even in the case of an object with a well 
resolved HMXB emission. The picture is emerging, the details 
are stiU missing. 
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APPENDIX 

In this Appendix we describe the reduction and analysis of the 
ASCA and BeppoSAX data for the nearby starburst galaxy M 82. 

— Data Reduction — 

ASCA data. The available archival spectral data for the galaxies M 82, 
M 83, NGC 253, NGC 2146, NGC 2903, NGC 3256, NGC 3310) were 
retrieved from the HEASARC for analysis. Net spectra for entire ob- 
servations, typically of duration 20 — 30 ks, were obtained, along with 
response matrices. Provided by the project were also the background 
spectra of the GIS detectors. The SIS backgrounds have been com- 
puted by accumulating cleaned counts arriving in a wide annulus sur- 
rounding (but excluding) the source into background spectra. Using 
software provided by the project we then generated response matrices 
for these accumulations Spectral fitting was performed simultane- 
ously on the two SIS and two GIS detectors, with a bandpass of 0.4-8 

Note that the derived spectral parameters are sensitive to the back- 
ground subtraction. Without a proper background subtraction, the hard 
X-ray spectrum of the source will appear systematically flatter, and the 
HMXB contribution will be overestimated. 



keV for SIS and 0.8-8 keV for GIS. We made sure we could closely re- 
produce the data and the results presented for the galaxies in the list by 
the original authors (M 82, NGC 253: Ptak et al. 1997; M 83: Okada 
et al. 1997; NGC 2146: Delia Ceca et al. 1999; NGC 2903: Mizuno 
et al. 1998; NGC 3256: Moran et al. 1999; NGC 3310: Zezas et al. 
1998). 

BeppoSAX data. The reduction of the BeppoSAX data for M 82 is de- 
scribed in detail in Cappi et al. (1999). Here let it suffice to say that 
our analysis is restricted to the 0.1 — 4.5 keV and 1.5 — 10 keV energy 
bands of the LEGS and MFCS instruments, where the response ma- 
trices are best calibrated. The extraction region, 4', has been chosen 
in order to maximize the S/N ratio. Standard blank-sky files (provided 
by the BeppoSAX Data Center) were used for background subtraction: 
the background turned out to provide ~ 7% of the total counts at 6 
keV. The detection of M 82 in the PDS data (nominal band: 13 - 300 
keV) is statistically significant up to ~30 keV. However, some degree 
of contamination from M 81 is possible due to the latter's higher flux 
(by a factor ~2.5) and angular proximity (which ensures that about 
half of its flux is collected). On the other hand, the spectral slopes of 
M 81 and M 82 in the PDS band do differ from each other, suggesting 
a tighter constraint on the contamination from M 81. To add uncer- 
tainty to the issue, it should be recalled that M 81 is known to have 
varied by a factor 4 on a timescale of months/years. These consid- 
erations led Cappi et al. (1999) to a conservatively assume that the 
PDS data provided only upper limits. In our analysis here, although 
we used all (including PDS) data, we made sure that the main result 
of our analysis (i.e., the estimate of the fractional 2-10 keV luminos- 
ity arising from the HMXB population) is substantially solid against 
including or rejecting the PDS data. 

— Data Analysis — 

Following Persic & Rephaeli (2002), the model spectrum we have 
used to analyze the M 82 emission comprises a thermal plasma plus 
the stellar-endproduct populations of SNRs, HMXBs, and LMXBs. 

As for the thermal plasma component, we found that the details 
of the fit at €^ 1 keV do not appreciably affect the results of the fit 
(i.e., the amplitudes of the stellar endproducts components) at €^ 2 
keV, and vice versa. However, an accurate low-energy fit was required 
in order to give an acceptable value of xt- So we used a multi-phase 
thermal plasma to fit the < 1 keV data, allowing both temperature and 
amplitude to vary, but forcing one same chemical abundance for all 
the phases. In fact, we do have theoretical insights (e.g., Suchkov et 
al. 1994) and some observational evidence (e.g., Dahlem et al. 1998, 
Strickland et al. 2000) on the existence of multi-phase X-ray emitting 
thermal plasmas in starburst environments, but no well-constrained 
ranges for the temperature and the chemical composition are available 
other than the rather generic (and expected, based on energy consider- 
ations) hint that the temperatures should be kT'^ 1 keV, and the chem- 
ical abundances should be Z^-, Zq. Specifically, we found that three 
plasma phases, described by Raymond-Smith (1977) models having 
temperatures kT = 0.065, 0.45, 0.75 keV and chemical abundance 
Z — 0.1 Zq, provide an acceptable fit to the low-energy data. The de- 
rived temperatures are within the range yielded for the X-ray emitting 
plasma by numerical simulations of the interaction between SN-driven 
winds and the ambient ISM in starbursts (Stickland & Stevens 2000). 
These parameter values were determined from an exploration of the 
parameter space < kTi < kT2 < fcTs < 1, with the constraint 
that Zi = Z2 = Z3. Once these parameters were determined, they 
were kept frozen (in the best-fit procedure) while the corresponding 
amplitudes Ai, A2 and A3 were left free to vary. As concerns the stel- 
lar end-products, we assumed the average spectral shapes observed 
for the corresponding Galactic populations to hold also in the case of 
M 82. 
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Table 4. Spectral fitting results for M 82(") . 
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The quoted errors represent the estimated 1 a confidence intervals, and are calculated by the X-ray spectral-fitting program XSPEC from 
the derivatives of the fit statistic with respect to the model parameters. 

Raymond-Smith (1977) plasma models with kTi = 0.065 keV, fcTa = 0.45 keV, fcTg = 0.75 keV, and Zi = Z2 = Z3 = 0.1 Zq. The 
amplitude of the RS plasma model is defined as A = 10"^*/47r((l + z) Da) J ngTiH, where Da is the angular size distance to the source 
(expressed in cm), rig and nn are the electron and hydrogen densities (in cm^^). 

HI column density, in excess of the foreground Galactic value, leading to photoelectric absorption of emission: units are 10^^ cm~^. 

RS plasma model with fcT = 2 and Z = Zq. 
^"^^ Power-law model of the form A{e) = A5{e/1 keV)~'^ where As is the amplitude at 1 keV (expressed in photons keV~^ cm~^ s"'^) and 
r = 1.2 is the photon index. 

Cutoff power-law model of the form A{€) = Aq (e/1 keV)~'^e~^/^°, where Ag is the amplitude at 1 keV (expressed in photons keV~^ 
cm~^ s~^), a = 1.4 is the photon index, and ec = 7.5 keV is the cutoff energy of the exponential cutoff (in keV). 
All model spectra are absorbed through a foreground Galactic HI column density Nm = 0.427 x 10^^ cm~^. 

The results of the fitting, for the separate as well as joint ASCA 
and BeppoSAX data sets, are shown in Table 4. Within the errors, the 
results of the separate fits and of the joint fit are mutually consistent. 



